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Abstract Electromagnetic fluctuation-induced forces between atoms and surfaces 
are generaUy known as Casimir-Polder interactions. The exact knowledge of these 
forces is rapidly becoming important in modern experimental set-ups and for tech- 
nological applications. Recent theoretical and experimental investigations have 
shown that such an interaction is tunable in strength and sign, opening new perspec- 
tives to investigate aspects of quantum field theory and condensed-matter physics. In 
this Chapter we review the theory of fluctuation-induced interactions between atoms 
and a surface, paying particular attention to the physical characterization of the sys- 
tem. We also survey some recent developments concerning the role of temperature, 
situations out of thermal equilibrium, and measurements involving ultra-cold atoms. 



1 Introduction 



In the last decade remarkable progress in trapping and manipulating atoms has 
opened a wide horizon to new and challenging experimental set-ups. Precision tests 
of both quantum mechanics and quantum electrodynamics have become possible 
through the capacity of addressing single trapped particles |[T][2| and of cooling ul- 
tracold gases down to Bose-Einstein condensation fS^Sl]. This stunning progress is 
also very profitable to other fundamental areas of physics and to technology. For 
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example, ultracold gases have been suggested as probes in interesting experimental 
proposals aiming at very accurate tests of the gravity law |[6}|8j, looking for ex- 
tra forces predicted by different grand-unified theories f9l (see also the chapter by 
Milton in this volume for detailed discussions on the interplay between Casimir en- 
ergy and Gravity). Technologically speaking, one paradigmatic example of this new 
frontier is provided by atom chips |10, 11|. In these tiny devices, a cloud of atoms 
(typically alkalis like Sodium, Rubidium or Cesium) is magnetically or optically 
trapped above a patterned surface, reaching relatively short distances between a 



few microns to hundreds of microns 1 12 - 14|. The micro-machined surface patterns 



form a system of conducting wires, which are used to control the atomic cloud by 
tuning an external induced current (also superconducting wires have been demon- 
strated IBHiTi). 

At a fundamental level, all these systems have in common to be strongly influ- 
enced by all kinds of atom- surface interactions. A particular category are fluctuation- 
induced forces, of which the most prominent representative is the van der Waals 
interaction 1 18|. These forces usually derive from a potential with a characteristic 
power-law dependence 

van der Waals limit : V = Vydw ^ ^ ' 

where L is the distance between the objects (two atoms or a surface and a atom) and 
the exponent depends upon physical parameters and geometry of the system (n = 3 
for an atom and a thick plate, see Fig[T]). Historically speaking, the existence of 
this kind of interaction was postulated long before it was experimentally possible to 
address single atoms [20J . The first quantum-mechanical theory was formulated by 
F. London in the thirties using the idea that the quantum mechanical uncertainty of 
electrons in atoms can be translated into fluctuating electric dipole moments \2T\ . 
London found that two atoms attract each other following ([T]) with an exponent 
n = 6. London's theory was extensively applied in studying colloidal suspensions 
|[22 | which provided confirmations of its validity but also showed its limitations. 

The next step was taken by H.B.G. Casimir and his student D. Polder [23 | who 
applied the framework of quantum electrodynamics, including the concept of vac- 
uum (field) fluctuations. They generalized the London- van der Waals formula by 
relaxing the electrostatic approximation, in other words, including the effect of re- 
tardation. The main success of Casimir-Polder theory was to provide an explanation 
for the change in the power law exponent observed in some experiments 1 22 1 . Indeed 
for distances larger than a characteristic length scale Ao of the system, the effect of 
retardation can no longer be neglected, and this leads to 

Casimir-Polder limit L > Aq: V = Vcp oc , (2) 

and therefore to the dependence typical of the Casimir-Polder interaction be- 
tween two atoms. The scale Aq is in this case the wavelength of the main atomic 
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Fig. 1 Atom-surface potential (free energy of interaction V{d)) vs. distance d between a ^^Rb atom 
and a SiOi (sapphire) substrate, multiplied by . The potential is calculated using the theory of 
Dzyaloshinskii, Lifshitz, and Pitaevskii 1 19 Ch.VIII]. Note the logarithmic scale and the sign. The 
figure, adapted from Fig. 3 of |4 |, shows the potential slIT = 300 K (solid line), at 7 = OK (dash- 
dotted line), and the three asymptotic behaviors (dotted lines): van der Waals-London oc —l/d^, 
Casimir-Polder (oc -l/d"^), and Lifshitz (oc -T/d^).) 



absorption lines, which is in the visible to near infrared for typical alkali atoms, a 
few hundreds of nanometers. 

The estimates ([lj|2]) apply at T = when only quantum fluctuations play a role. 
If the temperature is nonzero, another length scale comes into play, the thermal or 
Wien wavelength 

Ar = ^, (3) 

which corresponds to the wave length where the thermal radiation spectrum peaks. 
Calculations of the atom-surface interaction using thermal quantum field theory 
have been pioneered by Dzyaloshinskii, Lifshitz, and Pitaevskii |19,24|. They were 
able to recover the van der Waals and the Casimir-Polder potentials as limit behav- 
ior of a more general expression and to confirm, quite surprisingly, that at distances 
L <C the interactions are typically dominated by quantum fluctuations, the main 
reason being that their spectrum is much wider than that of the thermal field (which 
is constrained by the Bose-Einstein distribution jSsj). At distances Xt they 
show that the potential shows again a cross-over from the Casimir-Polder to the 
Lifshitz asymptote: 



Lifshitz limit L > Ar : V = Vl ^ 



vdW 



(4) 
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where (Oq is the (angular) frequency corresponding to Aq. This potential that scales 
with temperature is actually a free energy of interaction and is also known as the 
Keesom potential between polar molecules: there, the dipoles are rotating freely un- 
der the influence of thermal fluctuations 1 26] . In this case (Rydberg atoms provide 
another example), the particle resonances overlap with the thermal spectrum, and the 
Casimir-Polder regime is actually absent. Eq.(|4]) predicts an apparent enhancement, 
at nonzero temperature, of the fluctuation-induced interaction. This does not neces- 
sarily happen, however, because the molecular polarizabilities are also temperature- 
dependent |26}j28|. 

In Fig[T] above, we considered the case of an alkali atom whose peak absorption 
wavelength Aq is much shorter than the Wien wavelength. The Lifshitz tail is then 
much smaller than the van der Waals potential. Note that is of the order of a 
few micrometers at room temperature, comparable to the smallest atom-surface dis- 
tances achieved so far in atom chips. The crossover between the Casimir-Polder and 
the Lifshitz regimes can thus be explored in these set-ups. We discuss corresponding 
experiments in Sec. [5] 

In the following sections, we start with a derivation of the interaction between an 
atom and a general electromagnetic environment (Sec|2]). We will refer to it using the 
term "atom- surface interaction" or "Casimir-Polder interaction". This latter is also 
of common use in the literature to stress the fluctuation-induced nature of the inter- 
action, although the term "Casimir-Polder" more correctly indicates the potential in 
the retarded limit (see Eq.([2])). The result will be valid within a second-order per- 
turbation theory and can be easily adapted to specific geometries. We provide some 
details on a planar surface (Sec|3]). Situations out of global thermal equilibrium are 
discussed in Sec|4j dealing with forces on ultracold atoms in a general radiation en- 
vironment (the temperature of the surface and that of the surrounding environment 
are not necessarily the same), and with radiative friction. The final Sec |5] sketches 
experiments with atomic beams and ultracold samples. 



2 Understanding atom-surface interactions 

The interaction between atoms and between atoms and surfaces plays a fundamen- 
tal role in many fields of physics, chemistry and technology (see also the chap- 
ter of DeKieviet et al in this volume for detailed discussions on modern experi- 
ments on atom- surface Casimir physics). From a quantum-mechanical point of view, 
fluctuation-induced forces are not surprising and almost a natural consequence of 
the initial assumptions. Indeed the existence of fluctuations, even at zero tempera- 
ture, is one of the most remarkable predictions of the theory. Each observable corre- 
sponding to a physically measurable quantity can be zero on average but its variance 
will always be nonzero if the system is not in one of its eigenstates. When two quan- 
tum systems interact, the dynamics of the fluctuations becomes richer: each subsys- 
tem experiences, in addition to its own fluctuations, an external, fluctuating force. 
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This becomes particularly clear in the case of a po- 
larizable particle (atom or nanosphere) interacting 
with the vacuum electromagnetic field. In vacuum, 
the electromagnetic field fluctuates not only by it- 
self, but also because there are fluctuating sources 
for it, like the electric dipole moment of the particle. 
At the same time, the particle's dipole is not only 
fluctuating on its own, but is also responding to the 
fluctuations of the electromagnetic field ||29j|30j. As 
a result, when two atoms are brought nearby, they 
interact through their fluctuations mediated by the 
electromagnetic field. Similarly when a particle is in 
proximity of a macroscopic object, electric currents 
fluctuating inside the object and the fluctuations of 

the particle lead to a distance-dependent force. This second case is complicated by 
the fact of dealing with a macroscopic object and its quantum-mechanical descrip- 
tion. However, if the medium responds linearly to an electromagnetic perturbation, 
the fluctuation-dissipation theorem [31 1 provides a connection between the field's 
autocorrelation function and its macroscopic response (or Green function). 

We will use the previous considerations as a starting point for the derivation of 
the Casimir-Polder interaction between a surface and an atom or also a nanoparticle. 
We will follow Refs |27,,32}|34l. Although this is not the unique approach [23^,24) 



Fig. 2 A schematic representa- 
tion of the atom in the trap near 
a surface. 



35 -4^, it provides a physically transparent way to reach our final result. 



2,1 Energy of a polarizable particle in an electromagnetic field 

When a polarizable particle is introduced in an electric field, the change in energy 
of the system is given by ||43J 

{d(f)-E(ro,f)) ^ 

where, since we are working in the Heisenberg representation, all the operators are 
time dependent. From the thermodynamic point of view the previous quantity is a 
free energy and gives the amount of work that can be extracted from the system by 
moving the particle: in our (thermodynamic) convention a negative free energy will 
correspond to an attractive interaction (binding energy). 

The expectation value (• • •) is taken over the (initial) state of the non-coupled 
system; d is the (electric or magnetic) dipole operator and E the corresponding field 
operator, evaluated at the dipole position fq. We are implicitly assuming that the size 
of the particle is small enough to locally probe the electromagnetic field. The factor 
1 /2 in ([5]) arises from the fact that we are considering the energy of a linear polar- 
izable particle in an external field, rather than a permanent dipole | [43| . Note that 
the choice of a particular ordering does not seem to be necessary at this stage since 
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the dipole operator and the electric field operator commute. The symmetric order 
proves, however, to be particularly useful if we want to attach a physical meaning 
to each single contribution to the energy |29,30l, see Eq.([T5 ). 



The Hamiltonian of the coupled system can be in general written as // = //q + V , 
where Hq is the sum of the Hamiltonians of the two isolated subsystems and V de- 
scribes the interaction between them. Starting from this, the equation of motion for 
an operator A (^) can be written in the following integral form (Heisenberg picture) 

A{t) =A^'^\t)^l f&xe'^^^^'-'\yMx)] ^i^o(^-T) ^ 
n Jo 

where ^ = was chosen as initial time and the superscript free indicates that the 
operator evolves with respect to the Hamiltonian of the uncoupled system (//q), i.e. 

A^'''{t) = e-^^^'Ae-r^^^' . (7) 

Now, within first-order perturbation theory, eq.([6]) can be solved by replacing the 
operator A{t^) under the integral by its corresponding free evolution. If the inter- 
action Hamiltonian is bilinear like in the case of the electric dipole interaction, 
V = — d • E(ro), we get the following approximate expression for the dipole: 



'\t) = |_^dTQ[dfr-(0,dfr-(r-T)]e(T)') •Efr-(ro,r-T) , 



(8b) 



and similarly for the field: 



E(r,0«Efr-(r,0+E'"(ro,f), (9a) 
E'"(ro,f) = fjit (^^[Efr-(r,O,Efr-(ro,f-T)]0(T)^ ■d'''%t-x) . (9b) 

The T-integral runs effectively over T > (note the step function 0(t)) because in 
Eq.([6]), only times T > after the initial time are relevant (causality). In addition, we 
have set the upper limit to T = oo assuming that there exists a transient time after 
which the system behaviour becomes stationary. This time can be estimated from 



the system operators in Eqs.( |8b|[9b| ): the commutators are either c-number functions 
that die out for time arguments that differ by more than T^, or taking the expectation 
value, one gets subsystem correlation functions with as correlation time. 

We have therefore that, within the first order perturbation theory, the dipole in 
addition to its unperturbed evolution (d^^^^(^)) "responds" linearly (d^^(^)) to an ex- 
ternal perturbation (in this case the electromagnetic field). The same also happens 
to the electromagnetic field where now the dipole is the external source of pertur- 



bation. The term in parenthesis under the integrals (8b 9b), when evaluated over 
a particular state, is called susceptibility and contains the detailed physical infor- 
mation about the linear response of the system to the perturbation ||29l|30l . In the 



particular case of a dipole the Fourier-transform of the susceptibility tensor is the 
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polarizability 



/CX3 1 



where we have taken the expectation value for a given quantum state | a). In the time 



domain (see Eq.(8b)), the atomic susceptibility links the hermitean dipole operator 
to a hermitean field operator; hence it must be a real function. The polarizability, 
being its Fourier transform, therefore satisfies 

[^^(a))]* = ^^(-a)*) . (11) 

In addition, because of causality, Eq.|To|) implies that ^^{co) must be analytical in 
the upper-half of the complex a)-plane. 

Similar conclusions hold for the electromagnetic field. If we assume that the 
dynamics of the field and the surrounding matter (other than the atom) can be com- 
pletely described in terms of bosonic operators 139^ .44. 45J, the result of the com- 
mutator in ( [9b| ) is a c-number and the susceptibility does not depend on the state of 
the radiation. The linearity of the Maxwell equations then ensures that the result of 
eq.([9]) for a point-like dipole is correct to all orders. A simple identification leads to 
the following expression: 

r ^[Efr-(r,O,Efr-(ro,O)]0(r)e"»'dr = W(r,ro,ft)) . (12) 

J-oo n 

where W is the electric field Green tensor, solution to Maxwell's equation 

-VrXVrX G(r,ro,a)) + -^e(r,(o)G(r,ro,(o) = ^5(r-ro) / , (13) 

^ y 

where / is the identity tensor and e(r, (o) is the local dielectric function (here 
supposed to be a scalar for simplicity) of the matter surrounding the dipole. In con- 
clusion, in frequency space the induced quantities can be described in terms of the 
retarded response functions |[43| 



d^^(co) = ^((o) •Efr^^(ro,a)) , (14) 
Ei^(r,a)) = ^(r,ro,a))-dfr^^(a)), 

where the frequency dependence and causality allow for a temporal delay. This is 
slightly schematic because the polarizability tensor is defined only when the average 
is taken. 

Expressions ^ and ^ formalize the considerations made at the beginning of 
this section: both the dipole moment d(^) and the field E{t) can be split into two 
parts, the (free) fluctuating part describes the free intrinsic fluctuation, while the 
induced part arises in perturbation theory from the dipole coupling | |25| . Eq.([5]) be- 
comes 
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{di"(f);Efr^^(ro,f)) (Ei"(ro,f);dfr^^(f)) _ ^^^^ 

We assumed a factorized initial state in which each free evolution operator is zero 
on average and where the correlations between the fluctuating parts are entirely 
encoded in the linear response functions, setting the correlation between the freely 
fluctuating components to zero. This assumption would break down at higher orders 
of perturbation theory. Note that while in Eq.([5]) the total dipole and field operators 
(Heisenberg picture) commute at equal times, this is no longer true for their 'in' and 
'free' constituents in Eq.(p3]). The choice of the symmetric order (indicated by the 
semicolon) allows one to see each term of the previous expression as the result of 
the quantum expectation value of a Hermitian operator and therefore to attach to it 
a physical meaning f29][30j|32l . The first term on the right hand side of ( p3] ) can be 
seen as the contribution to the Casimir-Polder energy coming from the fluctuations 
of the vacuum field; the second will be called the self-reaction term since it arises 
from the interaction of the dipole with the field generated by the dipole itself. 



2.2 Equilibrium fluctuations 

Consider now a configuration at global thermal equilibrium, i.e. when both the 
dipole and the field are in a thermal state at temperature T. In this case we can ap- 
ply the fluctuation-dissipation theorem | 31|. This milestone of the linear response 
theory connects the correlation of a generic observable of a system in thermal equi- 
librium at temperature T with the imaginary part of the linear susceptibility which 
characterizes the response to a weak perturbation. In our case the theorem holds 
separately for the dipole and the field and we have |46 -48| 



{Ef'\v,(o)Ef\vM))T = 2nh5{(0 + (o')coth(^^ym[Gij{r,r,(o)] ,(16) 
{4'%(0)df^i(0'))T = 27r^5(«) + a)')coth Malj{(i>)\ , (17) 

where the symbol (• • • )r define the quantum and the thermal average and, according 
to ([To]), ^^((o) defines the atomic polarizability operator evaluated at temperature 
T (see Eq.([25]) below). The coth function in Eqs.([T6 17) arises from the symmet- 
rically ordered average of bosonic annihilation and creation operators in a state of 
thermal equilibrium: 

{a"^ a ^ aa'^)T = I ^2N{co) = coth (^^^) , (18) 

\2kBT J 



where A/^ is the Bose-Einstein distribution. Note the asymptotic limits 
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hco jl T<^hco/kB 
coth— > S Ok T , (19) 

2^^^ llf T:^hco/kB 

in the "quantum" (low-temperature) and "classical" (high-temperature) limits. 

The expression given in ([16]) can be directly reconnected with the currents fluc- 
tuating inside the media surrounding the dipole. For these currents Rytov's the- 
ory [ 49 1 predicts a correlation similar to ([17]) where the role of the polarizability is 
now played by the dielectric function | |46}j48| (see Sec |4] below). This picture also 
lends itself to a natural generalization where the bodies are assumed to be in "local 



thermal equilibrium" (see Sec 4.4 ). 



Note that the field correlations are needed at the same position fq. The Green 
function is, however, divergent in this limit due to its free- space contribution 

^(r,rQ,ft)) = ^^o(r , ro,a)) + y(r,ro,co) . (20) 

free space scattered 

The corresponding part of the free energy provides the Lamb shift of the internal 
levels of the dipole immersed in the electromagnetic field 137\. This contribution 
is position-independent and does not contain any information about the interaction 
between the bodies and the dipole. Therefore it can be safely "hidden" in the (renor- 
malized) energy levels of the atom. The physical information about the interaction 
is indeed contained only in the scattered part of the Green function 1 3£| . If the body 
happens to be a plane surface, it follows from symmetry that the result can only 
depend on the dipole- surface distance L and we can set ^(ro;ro, (o) = 



Combining Eqs.([T5-17 ), we finally obtain that the free energy of a polarizable 



particle at nonzero temperature T has the following general form (Einstein summa- 
tion convention) 

^""4/ , (21) 

V ^ / 

We have used the reality condition ( [TT] ), implying that the imaginary part of both 
polarizability and Green tensors are odd in (O. Eq.([2T]) coincides with the expression 
of the atom-surface interaction derived by many authors | [23|[24l[27]|32}j42| . It is 
often expressed in an equivalent form using the analyticity of a^{co) and W {L^ co) 
in the upper half of the complex frequency plane. Performing a Wick rotation in the 
complex frequency plane yields the so-called Matsubara expansion p4l[50| 

^(L,r) = -ksTY^ ajj{i^,)<^jimi^n) , (22) 



n=0 



where the Matsubara frequencies = 2nnkBT /h arise from the poles of the hyper- 
bolic cotangent, and the prime in the sum indicates that the ^ = term comes with a 
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coefficient 1/2. Both W^(i(^) and ^(L,i(^) are real expressions for > because 
ofEq.([TO]). 

These considerations conclude our first general analysis of the Casimir-Polder 
interaction. In the following section we will analyze the particle response function 
appearing in the previous formulation, namely the atomic polarizability, and men- 
tion also the case of a nanoparticle. 



2.3 Polarizability tensor 

The previous results can be used for the interaction of a surface with atoms, 
molecules, particles or in general any (small) object that can be described with good 
approximation in terms of a electric-dipole polarizability tensor. Here we are going 
to review the polarizability of an atom and of a nanoparticle. 



2.3.1 Atoms 

The polarizability tensor is determined by the transition dipole matrix elements and 
the resonance frequencies. For an arbitrary atomic state \a) it can be written as 

a- (CO) - y ^ (23) 

where df^ is the matrix element between the states \b) and \a) of the / component 
of the electric dipole operator and (0^^ = {Eh — Ea)/h the corresponding transition 
frequency. The introduction of an infinitesimal imaginary part shifts the poles of the 
expression into the lower part of the complex frequency plane (^cota — iO^). which 
is mathematically equivalent to the causality requirement. The tensorial form of the 
previous expression allows to take into account a possible anisotropic response of 
the atom to an electric field. A simplification of the previous expression can be 
obtained averaging over the different levels and directions so that the polarizability 
tensor becomes a^j = Sfja^^^ with the scalar function 

The polarizability is exactly isotropic when several excited sublevels that are degen- 
erate in energy are summed over, like the np^j^^ orbitals of the hydrogen-like series. 
When the atom is in thermal equilibrium, we have to sum the polarizability over the 
states \a) with a Boltzmann weight: 

-Ea/kBT 
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where Z is the partition function. In the limit T ^ 0, we recover the polarizability 
for a ground state atom. For a single pair of levels \a) and \h), this leads to the 
following relation between the state-specific and the thermalized polarizabities: 

a5(co)«a«.(co)tanh^. (26) 

This is mainly meant to illustrate the temperature dependence, otherwise it is a quite 
crude approximation. The reason is that the coupling to other levels makes the polar- 
izabilities afj and a^j differ quite a lot. Electronically excited states are much more 
polarizable due to their larger electron orbitals. 



2.3.2 Nanospheres 

Let us consider now the case where the atom is replaced by a nanosphere r34"5r521 . 
Indeed, if the sphere radius R is smaller than the penetration depth and the radiation 
wavelength, we can neglect higher order multipoles in the Mie expansion jSSj and 
consider only the electric and magnetic dipole (the sphere is globally neutral). 

In this long- wavelength limit, the Clausius-Mossotti relation ||43j|54j| provides 
the electric polarizability 

asph(co)=4;reo^3|^|^. (27) 

where £((0) is the (scalar) dielectric function of the sphere material. The nanosphere 
has also a magnetic polarizability that arises because a time-dependent magnetic 
field induces circulating currents (Foucault currents) |43|. This leads to a diamag- 
netic response | [55| 

/3sph(«) = ' - \\R' ■ (28) 

Both polarizabilities are scalars. For a metallic sphere, the electric polarizability 
goes to a positive constant at zero frequency, while the magnetic one vanishes there 
and has a negative real part at low frequencies (diamagnetism). 

For a qualitative comparison to an atom, one can estimate the oscillator strength 
| [25| , defined by the integral over the imaginary part of the polarizability. For the 
atom we have 

rdcolm aat((0) - and Tdcolm ^^.{(d) - ^ , (29) 

Jo n Jo n 

where the Bohr radius gq and the Bohr magneton jIb give the overall scaling of the 
transition dipole moments. The following dimensionless ratio allows a comparison 
between the two: 
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where (Xf^ = / {Ane^fic) ^ 1/137 is the fine structure constant. The electric oscil- 
lator strength clearly dominates in the atom. 

Let us compare to a metallic nanosphere (gold is often used in experiments) 
and assume a Drude model ( [ST] ) for the dielectric function. In terms of the volume 
V = 47lR^/3, we get an electric oscillator strength 



f 

Jo 



3 CO y 
I Imasph((o)J(0= -£b^V + ^(^). (31) 
2 V3 COp 



where COp/ \f3 is the resonance frequency of the particle plasmon mode (the pole of 
asph((o), Eq.([27|)). This is much larger than for an atom if the nanoparticle radius 
satisfies ao <C 7? <C A^, i.e., a few nanometers. The magnetic oscillator strength can 
be estimated as 

f .mft„.(«W<»=|^rlog(5=)(^)V. (32) 

where we took aOp as a cutoff frequency to make the integral convergent (at higher 



frequencies, Eq.( 28 ) does not apply any more). We have used the plasma wavelength 
Xp = Inc/cOp (~ 100 nm for gold). Similar to an atom, the nanoparticle response is 
dominantly electric, but the ratio of oscillator strengths can be tuned via the material 
parameters and the sphere size. The magnetic contribution to the particle- surface in- 
teraction is interesting because it features a quite different temperature dependence, 
see Ref. W\. 



2.4 Non-perturbative level shift 

In the previous section we saw that the main ingredient to derive the Casimir-Polder 
interaction between a particle and an object is the ability to solve for the dynamics 
of the joint system particle+electromagnetic field. Previously we limited ourself to 
a solution at the first order in the perturbation, implicitly motivated by the difficulty 
to solve exactly the dynamics of a multi-level atomic system coupled to a contin- 
uum of bosonic degrees of freedom (e.m. field). Things are different if we consider 
the linear coupling between two bosonic systems, i.e. if we describe the particle as a 
quantum harmonic oscillator. The linearity of the coupled system allows for an exact 
solution of its dynamics and even if the harmonic oscillator may be in some cases 
only a poor description of an atom [ 56 1, it is a good representation of a nanoparticle 
(the resonance frequency being the particle plasmon frequency). Generally, this ap- 
proach gives a first qualitative indication for the physics of the interaction p7}|6T| . 

The main idea we follow in this section is based upon a generalization of the 
"remarkable formula" of Ford, Lewis and O'Connell |62j|63J (see also |57}{6T|). 
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According to this formula, the free energy of a one-dimensional oscillator immersed 
in black body radiation is 

^FL0c(7^) = - r dco f{co,T)lm[dcolna{co)] , (33) 

71 Jo 

where /((O, T) is the free energy per mode, 



f{co,T)=kBTlog 



2sinh 



hco 

2kBf 



(34) 



and a{co) is the (generalized) susceptibility of the oscillator derived from Eq.(39) 
below. More precisely, ^{T) gives the difference between two free energies: the 
oscillator coupled to the radiation field and in equlibrium with it, on the one hand. 



and solely the radiation field, on the other. Eg. ([33]) is "remarkable" because the only 
system-relevant information needed here is the susceptibility function. 
In three dimensions, the polarizability becomes a tensor 

d((o) = ^(co) •£((()) , (35) 

where E{co) is the external electric field. In the case considered by Ford and 
O'Connell, there was no need to include a spatial dependence because of the ho- 
mogeneity and isotropy of the black body field. We are going to consider this sym- 
metry to be broken by the presence of some scattering object. As a consequence, the 
generalized susceptibility tensor becomes position-dependent W(a),ro). The spa- 
tial dependence is connected with the scattered part of the Green function and leads 
both to a position-dependent frequency renormalization and a damping rate. 

In order to get the expression of a(a), ro), let us consider for simplicity the equa- 
tion of motion of an isotropic oscillator with charge q interacting with the e.m. field 
near some scattering body (that is described by a dielectric constant). In frequency 
space, the (nonrelativistic) dynamics of the oscillator is described by 

m [(0^d((o) + coQd{co)] = ^^E(ro, co) , (36) 

where we have neglected the coupling with the magnetic field (first order in d/c). 
For the field we have 

V X V X E(r, co) ^£(co, r)E(r, co) = icojioUr, co) , (37) 

where the source current is j(r, co) = —icod{co)5{r — ro). 

Now, the formally exact solution for the operator E can be given in term of the 
(electric) Green tensor: 

E(r, CO) = Ef^^^(r, co) + ^(r, ro, co) • d{co) , (38) 
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where the Green tensor is the solution of Eq.( 13 ) given above. The field E^'^^r, co) 
is the electromagnetic field we would have without the oscillator and it is connected 
with the intrinsic fluctuations of the polarization field, or equivalently, of the currents 
in the body. Physically Eq. ([38]) states that the total electromagnetic field is given 
by the field present near the scattering object plus the field generated by the dipole. 
Introducing Eq.([38]) in Eq.([36|) we get 



-m(a)^ - (O^)d(a)) - ^^W(ro; ro, (o) • d((o) =^^E^'^^(ro, CO) . (39) 

The Green function ^(r;ro, co) solves an electromagnetic scattering problem and 
therefore, it decomposes naturally into a free-space field (as if the source dipole 
were isolated in vacuum), and the field scattered by the body, This is at the 
basis of the splitting in Eq.(20) discussed above. The free-space part ^ o{r;r() ^co) 



is a scalar in the coincidence limit because of the isotropy of space: part of the 
divergence (Re [Go]) can be reabsorbed into mass renormalization, mcOg ^ mcOQ, 
and part (Im[Go]) gives rise to dissipation (damping rate y{co)). Therefore Eq.([39|) 
can be rewritten as 

2 \ 2 

-(0^-i7((o)(0 + d)^-— ^(ro;ro,a)) ) •d((o) = — E^'^^(ro,(o) . (40) 
m J m 

The free electromagnetic field plays here the role of an external force and therefore 
the generalized (or "dressed") polarizability tensor is given by 

^((0,L) = av(co) (^1 -av(co)^(ro;ro,co)^ , (41) 
where we have defined 

^2 _^ 

av((o) = ^ (-0)^-17(0)) + . (42) 

If Ford, Lewis and O'Connell's result is generalized to a three-dimensional oscil- 
lator, a trace operation appears before the logarithm in Eq. ( [33] ). Using the identity 
tr log ^ = log det one gets 



1 f°° 

^FLOc{T) = - / da)/((0,r)Im[(9a)lndet^(a),J)] 

71 Jo 

dco /((0,r)lm [dco lnay{co)] 



1 r 

71 Jo 

1 /"°° 

71 Jo 



dco /((0,r)lm ^ojlndet (^1 - av((o)^(ro;ro, o))^ 



(43) 



The first term is distance-independent and coincides with the free energy of an iso- 
lated oscillator in the electromagnetic vacuum. It can be interpreted as a free- space 
Lamb shift. The second part of Eq. ( [43] ) is distance-dependent and therefore gives 
rise to the Casimir-Polder interaction. With help of a partial integration, we finally 
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( dcocothf I Im Tin detfl-av;( co)V(L, CO) 
2% Jq \2kBT J L V 



(44) 



The previous result can be easily generalized to the case of an anisotropic oscillator 
by just replacing the vacuum polarizability with the respective tensor. 



The usual expression (21) for the Casimir Polder free energy is recovered by 
assuming a weak atom-field interaction. Expanding the logarithm to first order we 
get 



^ = Tr / dco cothl —^]lm\ay(co)W m CO) 

In Jo \2kBT J 1^^ ^ ^ ' \ 



(45) 



From a scattering point of view, this approximation is equivalent to neglecting the 
multiple reflections of the electromagnetic field between oscillator and surface. At 
short distance to the surface, these reflections become relevant; the next-order cor- 
rection to the van der Waals interaction arising from ( [44| ) is discussed in Sec. 3. 4 of 
chapter by DeKieviet et al in this volume. 

Note that although very similar, Eqs.([2T]) and ( [45] ) are not identical. Eq.([2T]), 
applied to an oscillator atom, would have featured the bare polarizability 

a{co) = , ^^^"^ , (46) 
^ ^ (o2-(a) + i0+)2 ' 

where the infinitesimal imaginary part iO^ ensures causality. Eq.([45]) involves, on 
the contrary, the renormalized or vacuum-dressed polarizability which is causal by 
default. In other words, it contains a summation over an infinite subclass of terms in 
the perturbation series. 

Finally, as a general remark and in connection with the scattering interpretation 
of dispersion forces (see the chapters by Lambrecht et al. and by Rahi et al. in 
this volume for detailed discussions on the calculation of the Casimir effect within 
the framework of the scattering theory), within the theory of two linearly coupled 
linear systems, the susceptibilities involved in the description of the equilibrium 
Casimir-Polder interaction are the isolated and dressed ones (isolated scatters). This 
means that, within a linear response theory, or equivalently up to the first order in 
the perturbation theory, the susceptibilities are not modified by the presence of the 
other scatters but only dressed by the electromagnetic field. In our case, this means 
that 7(a)) or d) in Eqs.(|44|) or ( [45] ) do not depend on fq. 



3 Atoms and a planar surface 

Let us consider for definiteness the Casimir-Polder potential near a planar surface, 
with a distance L between the atom and surface. The Green function is in this case 
explicitly known and is given in the following subsection. 



16 



Francesco Intravaia, Carsten Henkel, and Mauro Antezza 



3.1 Behaviour of the Green function 

We sketch here the qualitative behaviour of the electromagnetic Green function near 
a planar surface that can be calculated analytically. Let the atom (source dipole) be 
on the positive z-axis at a distance L from a medium that occupies the half-space 
below the xj-plane. 



3.1.1 Reflection coefficients and material response 

The electric Green tensor ^(r,ro, co) is needed for coincident positions r = ro; by 
symmetry it is diagonal and invariant under rotations in the xj-plane p6||46||64[|65| : 



CX3 

^(L,co) / kdkK 





r™{(0,k)^^r'^{co,k)][xx- 



K*2 



where £o is the vacuum permittivity, /: = |k| is the modulus of the in-plane wave 
vector, and xx, yy, zz are the cartesian dyadic products. We consider here a local 
and isotropic medium, excluding the regime of the anomalous skin effect 1 66] . The 
Fresnel formulae then give the following reflection coefficients in the TE- and TM- 
polarization (also known as s- and p-polarization) |[43|: 



r^>,^) = 44^, r™(a),^) = 44^, (48) 

where K, are the propagation constants in vacuum and in the medium, respec- 
tively: 

, Km = yk^-£{C0)ll{C0)-^ . (49) 

The square roots are defined so that Im /c, Im /c^ < and Re /c, Re Km > 0. In par- 
ticular K is either real or pure imaginary. The corresponding frequencies and wave 
vectors define two regions in the {co^k) plane p4|: Evanescent region CO < ck: the 
electromagnetic field propagates only parallel to the interface and decays exponen- 
tially (k* > 0) in the orthogonal direction. Propagating region CO > ck: the electro- 
magnetic field also propagates (Re k* = 0) in the orthogonal direction. Note that the 
magnetic Green tensor can be obtained from the electric one by swapping the 
reflection coefficients \61j : 

Mo 
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All information about the optical properties of the surface is encoded in the re- 
sponse functions £{co) and il{co). For the sake of simplicity, we focus in the follow- 
ing on a nonmagnetic, metallic medium (jl{co) = 1) and use the Drude model [ |43| : 



col 



£{co) = l- , ^ . , , (51) 

where COp is the plasma frequency (usually for metals in the UV regime). The dis- 
sipation rate 7 takes account of all dissipative phenomena (impurities, electron- 
phonon scattering, etc.) in the metal | [68| and generally 7/0)^ <C 1 (~ 10~^ for gold). 



3.1.2 Distance dependence of the Green tensor 



The Drude model includes Ohmic dissipation in a very characteristic way, through 
the parameter 7. This affects the physical length scales of the system (see Ref. | [69| 
for a review). In our case the relevant ones are the photon wavelength in vacuum Ao, 
and the skin depth in the medium 5a). While the first is simply given by 



271C 
CO 



(52) 



the second is defined in terms of the low frequency behavior of the dielectric func- 
tion 



1 

5(0 



CO 



lm^/£{co) 



(for (0<C 7) 



(53) 



where D = yc^/cOp is the diffusion coefficient for the magnetic field in a medium 



with Ohmic damping | [43| . The skin depth gives a measure of the penetration of 
the electromagnetic field in the medium (~ 0.79/im at 10 GHz for gold). If we 
have 5a) <C Ao), the dependence of the Green function on L is quite different in 
the following three domains: (i) the sub-skin-depth region, L <C 5a), (ii) the non- 
retarded region, 5(o <^X(o, (iii) the retarded region: X(o <C L. In zones (i) and 
(ii), retardation can be neglected (van-der-Waals zone), while in zone (iii), it leads 
to a different power law (Casimir-Polder zone) for the atom-surface interaction. 
In the three regimes, different approximations for the reflection coefficients that 



appear in (47 ) can be made. In the sub -skin- depth zone | [67| , we have l/5a) ^ 
1 / Aa) and 



r^^(a),y^) ^ [e{co)-\\ 



TM 



{CO,k) 



CO^ 



4c^k^ ' 



e{co) 



co^ 



£{C0)-\ 

e(a)) + l ' e{co)^\c^k^ _ 



(54) 



At intermediate distances in the non-retarded zone, the wave vector is 1/Aa) <C ^ <C 
l/5a), hence 
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TE/ - ^ . . 2 ck 



r'^'ico.k) ^ -1+i 



r™(a),fc)«l+i^:^^. (55) 

Finally, in the retarded zone we can consider ^ <C l/Ao, <C so that 

rTE(co,^)^-l + — 

r™(a),^)^l-^^. (56) 



Note that the first terms in Eqs.( [55|[56| ) correspond to a perfectly reflecting medium 
(formally, e ^ oo). 

The asymptotics of the Green tensor that correlate to these distance regimes are 
obtained by performing the ^-integration in Eq.([47]) with the above approximations 
for the reflection coefficients. The leading-order results are collected in Table [T] 
One notes that the zz-component is larger by a factor 2 compared to the xx- and 
jj-components. This difference between the normal and parallel dipoles can be un- 
derstood by the method of images 143]. 

The magnetic response for a normally conducting metal in the sub- skin-depth 
regime is purely imaginary and scales linearly with the frequency (0\ the reflected 
magnetic field is generated by induction. A significant response to low-frequency 
magnetic fields appears for superconductors because of the Meissner-Ochsenfeld 
effect |70|. In contrast, the electric response is strong for all conductors because 
surface charges screen the electric field efficiently. 

The imaginary part of the trace of the Green tensor determines the local mode 
density (per frequency) for the electric or magnetic fields |71 1. These can be com- 
pared directly after multiplying by £o (or 1 /Mo). respectively. As is discussed in 
Refs. |[69j[7r|, in the sub-skin-depth regime near a metallic surface, the field fluc- 



tuations are mainly of magnetic nature. This can be traced back to the efficient 
screening by surface charges connected with electric fields. Magnetic fields, how- 
ever, cross the surface much more easily as surface currents are absent (except for 
superconductors). This reveals, to the vacuum outside the metal, the thermally ex- 
cited currents within the bulk. 



Table 1 Magnetic and electric Green tensors at a planar surface. In this case the other elements 
have the asymptotes = ^xx, = ^^xx, and similarly for ^u. The off-diagonal elements 
vanish. The expressions are for metals where \£{co) \ ^ 1. 



Sub- skin depth Non-retarded 



Retarded 



1 



£{(0)) 



?>2%5lL 



IM) 
' 327rL3 



3271 £oL^ V £((0) 

Mo 



327rL3 



1-i- 



1- 
47rL 



AtiL 1 


"47rL" 




^0) 2 


. ^(0 _ 





A7l\L/X(o 
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3.2 Asymptotic power laws 



To begin with, particle and field are both at zero temperature. The Matsubara series 
in Eq.([22|) can be replaced by an integral over imaginary frequencies: 

^ = ~hl d^I«I/-(i^)%(^'i^) . (57) 

^ 

where we have used the fact that the Green tensor is diagonal. Alternatively one 
can get the previous result by taking the limit T ^ of Eq.([2T]) and performing 
a Wick rotation on the imaginary axis. One of the main advantages of this repre- 
sentation is that all functions in ( [57] ) are real. For electric dipole coupling, one has 
a| (i(^ ) , %i (L, i(^ ) > 0, and we can conclude that the Eq.( 57 ) is a binding energy and 
corresponds to an attractive force (see the chapter by Capasso et al in this volume 
for detailed discussion on repulsive fluctuation-induced forces in liquids). 

Along the imaginary axis, the Green tensor is dominated by an exponential 



e-2<§L/c^ see Eq.(47). This exponential suppresses large values of t, and the main 
contribution to the integral comes from the region <c/ (2L). If this value is smaller 
than the characteristic frequency 12^, say, of the atom or of the nanoparticle (the 
lowest transition in eq. ([23])), the polarizability can be approximated by its static 
value. Assuming an isotropic polarizability we get the Casimir-Polder asymptote 

ae=c/{2Qe)) 

which is the well known expression for the atom-surface Casimir-Polder interac- 
tion f23l . At short distance the polarizability limits the relevant frequency range to 
^ < Qe- Therefore for L <C we can replace Green tensor by its short distance 
approximation (see Table [T]) where it becomes independent of ^ . We recover then 
the van der Waals asymptote 

L«A,: ^,,^«--^-|-3|d§af,„(i0. (59) 



Similar expressions hold for the interaction due to a fluctuating magnetic dipole, the 
behaviour becoming more complicated when the distance becomes comparable to a 
characteristic skin depth (see Eq.([53]) and Ref. |27 |). 

If we write the Matsubara frequencies as = ^tchc/Xt (^ = 0, 1, 2, . . .), the tem- 
perature may be low enough so that the limit Xt ^ L holds. Then all Matsubara 
frequencies are relevant, and if they are dense enough (A^ ^ A^), the effect of tem- 
perature is negligible. The series in ( [22] ) is then well approximated by the integral in 
( [57] ). In the opposite (high-temperature) limit, one has <C L so that the exponen- 
tial behavior of the Green tensor limits the series in ( [22] ) to its first term recovering 
the Lifshitz asymptote 
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(60) 



We still have an attractive force. Note, however, that this attraction is mainly due 
to the classical part of the radiation, as the same result would be obtained with a 
polarizable object immersed into the thermal field. 



4 Beyond equilibrium 

4.1 Overview 

The theory presented so far has mainly considered atom, field, and surface to be 
in a state of global thermal equilibrium, characterized by the same temperature T. 
When one moves away from these conditions, the atom-surface interaction assumes 
novel features like metastable or unstable states, driven steady states with a nonzero 
energy flux etc. We review some of these aspects here, since they have also appeared 
in recent experiments (Sec|5|. On the theoretical side, there are a few controversial 
issues that are currently under investigation | [72}j74| . 



We start with atoms prepared in non-thermal states: ground or excited states that 
decay by emission or absorption of photons, and with atoms in motion where fric- 
tional forces appear. We then consider field-surface configurations out of global 
equilibrium like a surface surrounded by a vacuum chamber at different tempera- 
ture. 



4.2 Atoms in a given state and field in thermal equilibrium 

The generalization of the Casimir-Polder potential to an atom in a definite state 
\a) can be found, for example, in Wylie & Sipe |36|, Eqs.(4.3, 4.4). Now, the 
fluctuation-dissipation theorem for the dipole, Eq.([T7]), does not apply, but pertur- 
bation theory is still possible, with the result (summation over repeated indices /, j) 

^{Lj) = -kBTj^afj(l^nWji{L^^n)^^ , (61) 

^=0 b 

where is the state-dependent polarizability . The dipole matrix elements 

are written df^ = {a\di\b). The thermal occupation of photon modes (Bose-Einstein 
distribution) is 

Note the second term in Eq.([6T]) that is absent in thermal equilibrium. It involves 
the absorption and (stimulated) emission of photons on transitions a ^ b to other 
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quantum states, and the thermal occupation number N{(Oija) evaluated at the Bohr 
frequency fiCOija = Ei^—Ea. For this reason, it can be called resonant part. The first 
term that was also present in equilibrium now features the state-dependent polariz- 
ability tensor ^^{i^n)- This is the non-resonant part of the interaction. 

For the alkali atoms in their ground state \a) = the Bohr frequencies E^g are 
all positive (visible and near-infrared range) and much larger than typical laboratory 
temperatures (equivalent to the THz range), hence the thermal occupation numbers 
N{cOha) are negligibly small. By the same token, the ground- state polarizability is 
essentially the same as in thermal equilibrium (i^n) because the thermal oc- 
cupation of the excited states would come with an exponentially small Boltzmann 
weight. The atom-surface interaction is then indistinguishable from its global equi- 
librium form and dominated by the non-resonant part. 

With suitable laser fields, one can perform the spectroscopy of atom-surface in- 
teraction of excited states |a) = |^). It is also possible to prepare excited states 
by shining a resonant laser pulse on the atom. In front of a surface, the sec- 



ond term in Eq.([6T]) then plays a dominant role: the transition to the ground state 
where a real photon is emitted is accompanied by an energy shift proportional to 
Re[^(L, —COeg)] . This resonant contribution can be understood in terms of the radi- 
ation reaction of a classical dipole oscillator |[36]|76] : one would get the same result 
by asking for the frequency shift of an oscillating electric dipole in front of a surface 
- a simple interpretation in terms of an image dipole is possible at short distances 
(where the ^-dependence of the reflection coefficients ( [48] ) can be neglected). This 
term is essentially independent of temperature if the transition energy E^g is above 

A more familiar effect for the excited state is spontaneous decay, an example 
for a non-stationary situation one may encounter out of thermal equilibrium. We 
can interpret the resonant atom-surface interaction as the 'reactive counterpart' to 
this dissipative process. Indeed, the spontaneous decay rate is modified relative to 
its value in free space by the presence of the surface. This can also be calculated 
in classical terms, leading to a modification that involves the imaginary part of the 
Green tensor ^(L, (0^^). (The free-space contribution has a finite imaginary 
part.) In fact, both the decay rate and the interaction potential can be calculated from 
a complex self-energy of which Eq.(61 ) is the real part in the lowest non- vanishing 
order of perturbation theory. 

What happens if the Bohr frequencies hcOha become comparable toksTl This ap- 
plies, for example, to optically active vibrational transitions and to atoms in highly 
excited states (Rydberg atoms) where the energy levels are closely spaced. It is ob- 
vious from Eq.([6T]) that the resonant term is subject to cancellations among "up" 
(Eij > Ea) and "down" transitions {Ey <Ea) with nearly degenerate Bohr frequen- 
cies: the occupation numbers N{cOha) SindN{cOya) differ in sign, while Re[^(L, co)] 
is even in (0. To leading order in the high-temperature limit, the resonant term be- 
comes 
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I 

b>a 



(Oba 



I 

b'<a 



^ab' 



(63) 
(64) 



where the notation b > a and b < a means summing over states with energies Ei, 
above or below Ea. This is proportional to the anharmonicity of the atomic level 
spectrum around Ea. It vanishes exactly for a harmonic oscillator and reduces sig- 
nificantly the coefficient linear in temperature in weakly anharmonic regions of the 
atomic spectrum |[28| . 



43 Moving atoms 

An atom that moves in a radiation field can be subject to a frictional force, as pointed 
out by Einstein in his seminal 1917 paper on the blackbody spectrum |77|. This 
force originates from the aberration and the Doppler shift between the field the 
atom "sees" in its co-moving frame, and the "laboratory frame". (The latter frame 
is actually defined in terms of the thermal distribution function of the radiation field 
that is not Lorentz-invariant. Only the field's vacuum state in free space is Lorentz- 
invariant.) In addition, electric and magnetic fields mix under a Lorentz transforma- 
tion so that a moving electric dipole also carries a magnetic moment proportional 
to d X V where v is the (center-of-mass) velocity of the dipole (the Rontgen current 
discussed in Refs. ||7§^). 



4.3.1 Black body friction 

The free-space friction force f(v, T) is given by | ]82{[83| : 





f(v,r) = -v-f^^ Uco - , (65) 



where a{(o) is the polarizability of the atom (in its electronic ground state) and 
the approximation of slow motion (first order in v/c) has been made. For atomic 
transitions in the visible range, this force is exponentially suppressed by the Boltz- 



mann factor ~ ^-^(Oeg/kBT ^^^^ winning against the prefactor 1/r in Eq.(65 ). The 
physics behind this effect is the same as in Doppler cooling in two counterpropagat- 
ing laser beams: the friction arises from the frequency shift in the frame co-moving 
with the atom that breaks the efficiency of absorbing photons with counter- and co- 
propagating momenta. Einstein derived the Maxwell-Boltzmann distribution for the 
atomic velocities by balancing this radiative friction with the momentum recoil in 
randomly distributed directions as the absorbed photons are re-emitted, which leads 
to Brownian motion in velocity space ||77||84j| . Conversely, assuming thermal equi- 
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Hbrium and the vaHdity of the Einstein relation between momentum diffusion and 
friction, one can calculate the (linear) friction tensor ^ in f(v) = — from the 
correlation function of the force operator |[83][85| : 

'^=^ /dT(F(^ + T)F(0), (66) 

CX3 

where ¥{t) is the force operator in the Heisenberg picture, the operator product 



is synmietrized (as in Sec 2.1 ), and the average (• • • ) is taken at (global) thermal 
equilibrium. One recognizes in Eq.([66|) the zero-frequency component of the force 
correlation spectrum. 

The motion of atoms in the radiation field plays a key role for laser cooling of 
ultracold gases. Although a discussion of laser-induced forces is beyond the scope 
of this chapter, the basic principles can be illustrated by moving away from global 
equilibrium and assigning temperatures 7^, 7> to atom and field. An ultracold gas, 
immediately after switching off the lasers, would correspond to Ta in the nanoKelvin 
range, while 7> = 300 K is a good assumption for the fields in a non-cryogenic labo- 
ratory apparatus. Dedkov & Kyasov calculated the separate contributions from fluc- 
tuations of the atomic dipole and the field, respectively. We follow here Ref. [86]. 
Qualitatively speaking, the fluctuating dipole experiences a force when it emits a 
photon; this force is nonzero and depends on velocity, even after averaging over 
all emission directions, because the emission is isotropic only in the rest frame of 
the atom. The absorption of photons from the fluctuating field is accompanied by 
photon recoil, and here isotropy is broken because the Doppler shift brings cer- 
tain directions closer to the resonance frequency. (The same principle is behind the 
so-called Doppler cooling in two counterpropagating beams.) The sum of the two 
contributions takes the form (adapted from Eq.(29) of Ref. | [86[ ) 

Pi rd^k 

f{^) = -^^^J -^{k-f){corima{co^ ' ^^^^ 

co' = 7v(a) + k-v) , (68) 

where 7v = (1 — v^/c^)~^/^ is the relativistic Lorentz factor, and v, k are unit vectors 
along the atom's velocity and the photon momentum. The photon frequency in the 
"blackbody frame" (field temperature 7» is CO = c|k|, and N{co^T) is the Bose- 
Einstein distribution for a mode of energy quantum hco at temperature T. The term 
with N{co^Tf) gives the force due to absorption of thermal photons, while N{co\Ta) 
gives the force due to dipole fluctuations. The absorbed power has to be calculated in 
the atom's rest frame: the energy hco' times the photon number provides the electric 
field energy density in this frame, and the absorption spectrum Im [co^a{co^)] must be 
evaluated at the Doppler- shifted frequency co\ This shifted spectrum also appears 
in the fluctuation-dissipation theorem ^Tl) now applied locally in the atom's rest 
frame, and determines the dipole fluctuations. 
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At equilibrium and in the non-relativistic limit, the difference between the Bose- 
Einstein distributions can be expanded to give 

N{(0,Tf)-N{co',Ta) « -{k-y)do,N{(0,T) = (k-v)— ^J^^— — , (69) 

4smh (Aico/z^^r) 



and performing the angular integration, one recovers Eq.([65]). As another example, 
let us consider a ground-state atom (7a = 0) moving in a "hot" field 7> > with a 
small velocity. We can then put N{(o'^Ta) =0 since for free space photons, 0)' > 
(the positive-frequency part of the light cone in the ((0,k) space is a Lorentz- 
invariant set). Expanding in the Doppler shift, performing the angular integration 
and making a partial integration, one arrives at 

h 7 

f(v) ^ v ^^^^ /d(0(O^Ima^(a))(9a)[(«^A^(^<),7»] . (70) 

71 £oC ^ 

Note that the function co^N{(0,Tf) has a positive (negative) slope for CO < ksTp /h 
(CO > kBTf/h ), respectively. The velocity-dependent force thus accelerates the par- 
ticle, V • f > 0, if its absorption spectrum has a stronger weight at sub-thermal fre- 
quencies. This may happen for vibrational transitions in molecules and illustrates the 
unusual features that can happen in non-equilibrium situations. Drawing again the 
analogy to laser cooling, the radiative acceleration corresponds to the "anti-cooling" 
set-up where the laser beams have a frequency CO > COeg ("blue detuning"). Indeed, 
we have just found that the peak of the thermal spectrum occurs on the blue side of 
the atomic absorption lines. 



4.3.2 Radiative friction above a surface 

Near a surface, the fluctuations of the radiation field are distinct from free space, and 
are encoded in the surface-dependent Green function ^(L, (o), see Eqs.( 16 



20). In 



addition, one has to take into account that the available photon momenta differ, 
since also evanescent waves appear whose ^-vectors have components larger than 
co/c. All these properties can be expressed in terms of the electromagnetic Green 
tensor, assuming the field in thermal equilibrium. Let us consider for simplicity an 
atom with an isotropic polarizability tensor atj = a 5ij, moving at a non-relativistic 
velocity v. From Eq.(118) in Ref. | ,87J , one then gets a friction force 

f(v) = --4^ /d(o(-(9a)A^((0,r))Ima(a))(v-V)(vV')trIm^(r,r',a)) , (71) 





where the spatial derivatives are taken with respect to the two position variables of 
the Green tensor, and r = r' = rA{t) is taken afterwards. This expression neglects 
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terms of higher order in a that appear in the self-consistent polarizabiHty ( [4T] ); it 
can also be found from Eqs. (25,26) in Zurita Sanchez et al. | [83| . 

Note that friction is proportional to the local density of field states, encoded in 
the imaginary part of the (electric) Green tensor If the motion is parallel to a 
plane surface, the result only depends on the distance L and is independent of time. 
The friction force is comparable in magnitude to the free- space result ( [67] ) if the 
distance L is comparable or larger than the relevant wavelengths c/ft): the deriva- 
tives in Eg.!?!} are then of the order (v •¥)(¥• V) ~ |v| (co/c)^. At sub-wavelength 
distances, the non-retarded approximation for the Green tensor can be applied (see 
Table [!}, and the previous expression becomes of the order of |v|/L^. The remain- 
ing integral is then similar to the temperature-dependent part of the atom-surface 
interaction discussed in Sec l4.4.2l below. 



An expression that differs from Eq.(71) has been found by Scheel and Buh 



mann jSTJ who calculated the radiation force on a moving atom to first order in 
the velocity, and at zero temperature. Their analysis provides a splitting into reso- 
nant and non-resonant terms, similar to Eq.([6T]). For the ground state, the friction 
force is purely non-resonant and contains a contribution from the photonic mode 
density, similar to Eq.jfTT), and one from the Rontgen interaction that appears by 
evaluating the electric field in the frame co-moving with the atom. Another non- 
resonant friction force appears due to a velocity-dependent shift in the atomic reso- 
nance frequency, but it vanishes for ground- state atoms and for the motion parallel 
to a planar surface. The remaining friction force becomes in the non-retarded limit 
and for a Drude metal 

f(v) « ^ y l^^^l'^--^- (72) 

where the sum a > g is over excited states, the relevant dipole matrix elements are 
|^fl^|2 _ l^f^^p, Ta is the radiative width of the excited state (which also depends 
on L), and (0^ the surface plasmon resonance. The fluctuations of the electromagnetic 
field are calculated here without taking into account the "back-action" of the atom 
onto the medium (see Ref. 1 85 1 and the discussion below). 

We briefly mention that the behaviour of friction forces in the limit of zero tem- 
perature ("quantum friction") has been the subject of discussion that is still contin- 
uing (see also the chapter by Dalvit et al. in this volume for further discussion on 
quantum friction). An early result of Teodorovich on the friction force between two 
plates, linear in v with a nonzero coefficient as T ^ | [88| , has been challenged by 
Harris and Schaich [ 85 1. They point out that a charge or current fluctuation on one 
metallic plate can only dissipate by exciting electron-hole pairs in the other plate, 
but the cross-section for this process vanishes like . This argument does not hold, 
however, for Ohmic damping arising from impurity scattering. In addition, Ref. ["HFl 
points out that the fluctuations of the atomic dipole should be calculated with a po- 
larizabiHty that takes into account the presence of the surface. This self-consistent 
polarizabiHty has been discussed in Sec |2.4| and reduces the friction force, in par- 
ticular at short (non-retarded) distances. Carrying out the calculation for a metallic 
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surface and in the non-retarded regime, Harris & Schaich find the scahng 

r(v)-vi!|rJ5ffl^V, ,73) 



\47l£oCOs 

where afs is the fine structure constant, a(0) is the static polarizabiHty of the atom, 
e the elementary charge, and cOs = COp/V^ the surface plasmon frequency in the 
non-retarded limit. Note the different scaling with distance L compared to Eq.([72|). 



4.4 Nonequilibrium field 

A radiation field that is not in thermal equilibrium is a quite natural concept since un- 
der many circumstances, an observer is seeing radiation where the Poynting vector 
is nonzero (broken isotropy) and where the frequency spectrum is not given by the 
(observer's) temperature. The modelling of these fields can be done at various lev- 
els of accuracy: "radiative transfer" is a well-known example from astrophysics and 
from illumination engineering - this theory can be understood as a kinetic theory 
for a "photon gas". It is, in its simplest form, not a wave theory and therefore not 
applicable to the small length scales (micrometer and below) where atom-surface 
interactions are relevant. "Fluctuation electrodynamics" is a statistical description 
based on wave optics, developed by the school of S. M. Rytov 1 19, 89], and similar 
to optical coherence theory developed by E. Wolf and co-workers |25|. The main 
idea is that the radiation field is generated by sources whose spectrum is related to 
the local temperature and the material parameters of the radiating bodies. The field 
is calculated by solving the macroscopic Maxwell equations, where it is assumed 
that the matter response can be treated with linear response theory (medium permit- 
tivity or dielectric function e(x, (o) and permeability /i(x, (o)). This framework has 
been used to describe the quantized electromagnetic field, as discussed by Knoll and 
Welsch and their co-workers | [90| , by the group of Barnett | [9T| , see also the review 
paper Ref . [ 92 1 . Another application is radiative heat transfer and its enhancement 
between bodies that are closer than the thermal (Wien) wavelength, as reviewed in 
Refs. |[87][93j. The non-equilibrium heat flux between two bodies at different tem- 
peratures is naturally calculated from the expectation value of the Poynting vector. 

In this section we review the atom- surface interaction in the out-of-equlibrium 
configuration similar to the one studied in | [94| : the atom is close to a substrate hold 
at temperature 7^-, the whole being enclosed in a "cell" with walls (called "envi- 
ronment") at temperature Te. In the following we will only consider the electric 
atom-surface interaction and we will use the zero temperature atomic electric polar- 
izabiHty. In fact, the electric dipole transitions are mainly in the visible range and 
their equivalent in temperature (10^~^K) are not achieved in the experiment. There- 
fore the atom does not participate in the thermal exchange and can be considered in 
its ground state. 
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4.4.1 Fluctuation electrodynamics and radiative forces 

A very simple non-equilibrium situation occurs when an atom is located near a 
"heated body" whose temperature is larger than its "surroundings". As mentioned 




Fig. 3 Sketch of an atom- 
surface system with the field 
being out of thermal equilib- 
rium. Ts is the temperature of 
the substrate and Te is the tem- 
perature of the walls of the cell 
surrounding the atom- substrate 
system. If 7^ > Te, there is a 
nonzero radiative heat flux from 
the surface into the surrounding 
environment. 



above, it is quite obvious that the Poynting vector of the radiation field does not 
vanish: there is radiative heat flux from the body into the surrounding space. This 
flux is accompanied by a radiative force on the atom that depends on the atomic 
absorption spectrum, but also on the angular distribution of radiated and re-scattered 
photons. The atom-photon interaction, in this case, does not derive from the gradient 
of a potential. The basic concept is that of the radiation force F; it is given by (431 



F(r) = {diVEi{r)) + {l^iVBi{v)) + higher multipoles , (74) 

where we have written out only the contributions from the electric and magnetic 
dipole moments and the atom is assumed at rest at position r. (The generalization 



to a moving atom leads to the velocity-dependent forces discussed in Sec 4.3 ) As 



a general rule, the electric dipole is the dominant contribution for atoms whose 



absorption lines are in the visible range. Eq.(74) can be derived by averaging the 
Coulomb-Lorentz force over the charge and current distribution in the atom, assum- 
ing that the atomic size is small compared to the scale of variation (wavelength) of 
the electromagnetic field. The average (...) is taken with respect to the quantum 
state of atom and field, and operator products are taken in symmetrized form. 

The radiation force ( |74| ) can be evaluated with the scheme outlined in Sec |2.1| 
where the operators d and E(r) are split into "freely fluctuating" and "induced" 
parts. Carrying this through for the contribution of field fluctuations, leads to an 
expression of the form 

(J,-V£f-(r)) = ^ «,;(co)(4-(r, co) V£f-(r, co')) , (75) 
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where the spatial gradient of a field autocorrelation function appears. In a non- 



equilibrium situation, the fluctuation-dissipation theorem of Eq.( 16) cannot be ap- 
plied, and this field correlation must be calculated in a different way. In a similar 
way, one gets 

(Jf-V£™(r)) = / ^^S^,G,j{r,r,co'){df^^co)dj^%co')) , (76) 

where the gradient ^iGfj is evaluated with respect to the first position coordinate of 
the Green function. This term requires some regularization because of the divergent 
Green function at coincident positions. The correlation function of the atomic dipole 
can be calculated in its stationary state which could be a thermal equilibrium state 



or not, as discussed in Sec |4.2[ In the case of an ultracold atomic gas, it is clear 
that the atom can be at an effectively much lower temperature compared to the 
macroscopic bodies nearby. This is consistent with the perturbation theory behind 
the operator splitting into fluctuating and induced parts. In global equilibrium, when 
both fluctuation spectra are given by the fluctuation-dissipation theorem (16 T7|, it 



can be seen easily that the force reduces to the gradient of the equilibrium interaction 
potential ( [2T] ). 

Within Rytov's fluctuation electrodynamics, the fluctuating field is given in terms 



of its sources and the macroscopic Green function. Generalizing Eq.( 15 ), one gets 



^^^(r, (o) = y d^r Gij{r, r\ co)Pj{r\ co) + magnetization sources , (77) 

where we have not written down the contribution from the magnetization field 
that can be found in Ref. f95|. The polarization density Pj{r',co) describes the 
excitations of the material (dipole moment per unit volume). If the material is 
locally stationary, the polarization operator averages to zero and its correlations 
(/^(r, co)Pj{r\ co')) determine the field spectrum 

{Ei{r,co)Ej{r\co')) = J d^xdVG,-;^(r,x, a))Gj/(r',x', (oO(/\(x, (o)P/(x', co')) , 

(78) 

Making the key assumption of local thermal equilibrium at the temperature T{v), 
the source correlations are given by the local version of the fluctuation-dissipation 
theorem |[89j: 

(79) 

where e(r, O)) is the (dimensionless) dielectric function of the source medium, 
giving the polarization response to a local electric field. The assumption of a lo- 
cal and isotropic (scalar) dielectric function explains the occurrence of the terms 
Sij5{v — r'); this would not apply to ballistic semiconductors, for example, and to 
media with spatial dispersion, in general. The local temperature distribution r(r) 
should in that case be smooth on the length scale associated with spatial disper- 



Fluctuation-induced forces between atoms and surfaces: the Casimir-Polder interaction 29 

sion (Fermi wavelength, screening length, mean free path). Similar expressions for 
random sources are known as "quasi-homogeneous sources" and are studied in the 
theory of partially coherent fields | 25 , §5.2]. 
If we define a polarization spectrum by 

^p(r,(o) =;icoth(^^^^^^Im[eoe(r,co)] , (80) 

the field correlation function (TTSI) becomes 



(81) 

where Eq.(pT|) has been applied to the Green function. This expression was named 
"fluctuation-dissipation theorem of the second kind" by Eckhardt |47 | who ana- 
lyzed carefully its limits of applicability to non-equilibrium situations. Note that 



even if the sources are spatially decorrelated (different points r and r' in Eq.(79), 
the propagation of the field creates spatial coherence, similar to the lab class diffrac- 
tion experiment with a coherence slit. The spatial coherence properties of the field 
determine the order of magnitude of the gradient that is relevant for the radiation 



force in Eq.(75 ). 



Let us focus in the following on the correction to the atom- surface force due to 
the thermal radiation created by a "hot body". We assume that the atom is in its 



ground state and evaluate the dipole fluctuation spectrum in Eq.([76|) at an atomic 
temperature 7^ = 0. To identify the non-equilibrium part of the force, it is useful 
to split the field correlation spectrum in Eq.([75]) into its zero-temperature part and 
a thermal contribution, using coth{h CO /IksT) = sign((o)[l -^2N{\co\,T)] with the 
Bose-Einstein distribution N{co^T). The Rytov currents are constructed in such a 
way that at zero temperature, the fluctuation-dissipation theorem ([16]) for the field 
is satisfied. This can be achieved by allowing formally for a nonzero imaginary 
part Ime(r, ft)) > everywhere in space 1^48^20], or by combining the radiation 
of sources located inside a given body and located at infinity |l92l|96l. The two 



terms arising from dipole and field fluctuations then combine into a single expres- 



sion where one recognizes the gradient of the Casimir-Polder potential Eq.(21 ). This 



is discussed in detail in Refs. |33j[97|. The remaining part of the atom-surface force 



that depends on the body temperature is discussed now. 



4.4.2 Radiation force near a hot body 

Let us assume that the body has a homogeneous temperature T{x) = Ts and a spa- 
tially constant dielectric function. Using Eq.([8T]) and subtracting the T = limit, the 
spectrum of the nonequilibrium radiation (subscript 'neq') can then be expressed by 
the quantity 
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{Ei{r,co)Ej{r\(o'))r,,^ = 27i5{co ^ (0^)N{\colTs)nSij{ry ,co) , (82) 

Sij{ry,co)=lm[£o£{\co\)] j &\Gl(y,x,(o)Gjk{v' ,x,(o) , (83) 

s 

where the space integral is over the volume of the body. The tensor 5'/y(r,r', (o) 
captures the material composition of the body and its geometry relative to the ob- 
servation points. 



Referring to the force due to field fluctuations in Eq.(75 ), let us assume for sim- 
plicity that the atomic polarizability is isotropic, (Xij{co) = 5ija{co). We combine the 
integrand over positive and negative frequencies to isolate dispersive and absorptive 
contributions (ft) > 0) 

a*((o)V25';7(r,r,(o) + (o) -co) 

= Re[a(co)]Vr[^,7(r,r,(o)]+2Im[a((o)]Im[V2^,7(r,r,(o)] , (84) 

where V2 is the gradient with respect to the second coordinate of Sij, while Vr 
differentiates both coordinates. This form highlights that the non-equilibrium force 
separates in two |98j|99| contributions that are familiar in laser cooling |10Q||1QT| : 
a dipole force equal to the gradient of the electric energy density (proportional to 
5//(r,r, 0)) > 0). This is proportional to the real part of a and can be interpreted 
as the polarization energy of the atom in the thermal radiation field. The second 
term in Eq.(84) gives rise to radiation pressure, it is generall}^ proportional to the 



atomic absorption spectrum and the phase gradient of the field. The phase gradient 
can be identified with the local momentum of the emitted photons. By inspection 
of Eq.([83]) for a planar surface, one indeeds confirms that the force pushes the atom 
away from the thermal source. An illustration is given in Fig|4]for a nanoparticle 
above a surface, both made from semiconductor. The dielectric function e{(o) is 
of Lorentz-Drude form and uses parameters for SiC (see Ref. |33|). The arrows 
mark the resonance frequencies of transverse bulk phonon polaritons a>r and of the 
phonon polariton modes of surface {(0\) and particle {Ci>i). 

The radiation pressure force is quite difficult to observe with atomic transitions in 
the visible range because the peaks of the absorption spectrum are multiplied by the 
exponentially small Bose-Einstein factor N{cOag^Ts), even if the body temperature 
reaches the melting point. Alternative settings suggest polar molecules or Rydberg 
atoms |[28| with lower transition energies. In addition, some experiments are only 
sensitive to force gradients (see Sec|5]), and it can be shown that the radiation pres- 
sure above a planar surface (homogeneous temperature, infinite lateral size) does 
not change with distance. 

For the evaluation of the dipole force, the same argument related to the field 



temperature can be applied so that the atomic polarizability in Eq.(84) is well ap 



Note that when the dressed polarizability is used instead of the bare one, the polarizability has 
an imaginary part even in the absence of absorption (see discussion at the end of section [^J. 
This is equivalent to include the effects of the "radiative reaction" in the dynamic of the dipole 
(98]|102||1C6| as required by the conservation of energy and the optical theorem. 
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Fig. 4 (left) Spectrum of thermal radiation pressure exerted on a small spherical particle above a 
planar substrate (positive = repulsive, does not depend on distance). The arrows mark the substrate 
and particle resonances at Ree~^ (a>r) = (bulk phonon polariton) and Re8(ft)i 2) = — 1, — 2 (sur- 
face and particle phonon polariton). The force spectrum is given by hN{co,Ts) times the second 
term in Eq.|84), and normalized to (16 /3)hkL(kLa)^N((0,Ts) where a is the particle radius and 
kL = col/c the wavenumber of the longitudinal bulk polariton (Re8(ft)L) = 0). 
(right) Theoretical calculation of the atom- surface force, in and out of thermal equilibrium, taken 
from Ref. |94|, Fig. 2. The atom is Rubidium 87 in its electronic ground state, the surface is made 
of sapphire (SiOi). Note the strong variation of the non-equilibrium force, both in magnitude and 
in sign. A negative sign corresponds to an attractive interaction. 



proximated by its static value, Rea(()L)) ^ a(0). We thus find 
F^iP (r, 7i, = 0) = -VU^^^iv, Ts, Te = 0), 

U^^{r,Ts,TE=0) = -a(0) j^hN{coJs)Su{r,r,co) . 



(85) 
(86) 



4.4.3 General non equilibrium configuration and asymptoric behaviours 



In the general case both and Te can be different from zero |94|. The total force 



will be the sum of Fneq(r, 7$, Te = 0) given in the previous expression and of 

F^iP (r, Ts = ^Je)= F^jf (r, Te) - (r, 7> , = 0) , 



(87) 



that is the difference between thermal force at equilibrium at the temperature Te 
and the force in Eqs.([85]) and ([86]) with T^- and Te swapped. An illustration of the 
resulting force is given in Fig|4jright) for a planar surface at temperature Ts > Te. 
A large-distance asymptote of the non-equilibrium interaction can be derived in the 
form 1941 1 10411 1051 



32 



Francesco Intravaia, Carsten Henkel, and Mauro Antezza 



"-^^T^M^- ^neq(^,r.,r.) - -J^m ^^^^^ _ ^ , (88) 

where e(0) < oo is the static dielectric constant. The previous expression is vaUd 
for dielectric substrates, see Ref. | [94| for Drude metals. For Ts = Te, this formula 
vanishes, and one ends with the "global equilibrium" result of Eq.([6T]). 

Expression ( [88] ) shows that the configuration out of thermal equilibrium presents 
new features with respect to the equilibrium force. Indeed the force scales as the 
difference of the square of the temperatures and can be attractive or repulsive. For 
Te > Ts the force changes sign, going from attractive at small distance to repulsive 
at large distance (i.e. featuring a unstable equilibrium position in between), and it 
decays slower than the equilibrium configuration (oc L~^), leading therefore to a 
stronger force. This new feature was crucial for the first measurement of the thermal 
component of the surface-atom surface (see next section). Moreover, when a gas of 
atoms is placed in front of the surface, the non-equilibrium interaction can lead to 
interesting non-additive effects [97^^1Q6J . 



5 Measurements of the atom-surface force with cold atoms 
5.1 Overview 

We do not discuss here the regime of distances comparable to the atomic scale where 
atomic beams are diffracted and reveal the crystallography of the atomic structure 
of the surface. We consider also that the atoms are kept away even from being ph- 
ysisorbed in the van der Waals well (a few nanometers above the surface). One 
is then limited to distances above approximately one micron (otherwise the attrac- 
tive forces are difficult to balance by other means), but can take advantage of the 
techniques of laser cooling and micromanipulation and use even chemically very 
reactive atoms like the alkalis. 

The first attempts to measure atom-surface interactions in this context go back 
to the sixties, using atomic beam set-ups. In the last 20 years, technological im- 
provements have achieved the sensitivity required to detect with good accuracy and 
precision tiny forces. As examples, we mention the exquisite control over atomic 
beams provided by laser cooling |107||108 |, spin echo techniques that reveal the 
quantum reflection of metastable noble gas beams | (109[|110| (see also the chapter 
by DeKieviet et ah for detailed discussions on this topic), or the trapping of an ul- 
tracold laser-cooled gas in atom chip devices p] |llH . In this section we will briefly 
review some of the experiments which exploited cold atoms in order to investigate 
the Casimir-Polder force. 

By using different techniques, it has been possible to measure the atom-surface 
interaction (atomic level shift, potential, force, or force gradient, depending on the 
cases) both at small (0.1 /im < L < Aopt ~ 0.5/im) and large (L > 1/im) distances. 
Due to the rapid decrease of the interaction as the atom-surface separation becomes 



Fluctuation-induced forces between atoms and surfaces: the Casimir-Polder interaction 



33 



larger, the small-distance (van der Waals-London, Eq.([59|)) regime at L < Aopt is 
somewhat easier to detect. Recall that in this limit, only the vacuum fluctuation 
of electromagnetic field are relevant, and retardation can be ignored. More recent 
experiments explored the weaker interaction in the Casimir-Polder regime (58), 
Aopt <L < Xt, where retardation effects are relevant, but thermal fluctuations still 
negligible. Also the Lifshitz-Keesom regime at L > has been explored, where 
thermal fluctuations are dominant [see Eq.([60|)]. The theory of Dzyaloshinskii, 
Lifshitz, and Pitaevski (DLP, |19, Ch.VIII]) encompasses the three regimes with 
crossovers that are illustrated in Fig [T] for a Rubidium atom and a room temperature 
sapphire surface. 



5.2 From van der Waals to Casimir-Polder: equilibrium 

Typically, experiments have been performed at room temperature, and at thermal 
equilibrium, and used several techniques to measure the interaction, usually of me- 
chanical nature. 

The van der Waals-London regime has been first explored by its effect on the 
deflection of an atomic beam passing close to a substrate jll2f|116| . Such kind 
of experiments were almost qualitative, and hardly in agreement with the theory. 
Subsequently, more accurate measurements of the atom-surface interaction in this 
regime have been done by using dielectric surfaces "coated" with an evanescent 
laser wave that repels the atoms (atom mirrors, y_17j|), by atom diffraction from 
transmission grating s |118 119| , by quantum reflection | |1Q9||11Q| , and by spectro- 



scopic studies 1 107, 12Q| . 

The Casimir-Polder regime, where vacuum fluctuations of the electromagnetic 
fiel d an d the finite speed of light are relevant, was first studied experimentally 
in Here the force has been measured through an atomic beam deflection 

technique, which consists in letting an atomic beam (Na atoms in their ground state) 
pass across in a cavity made by two walls (gold plates), as one can see in fig|5] The 
atoms of the beam are drawn by the Casimir-Polder force to the walls, whose inten- 
sity depends on the atomic position within the cavity. Part of the atoms are deflected 
during their path and stick to the cavity walls without reaching the end of the cavity. 
Only few atoms will pass the whole cavity, and their flux is measured and related 
to the atom-surface interaction in the cavity. Such a measurement is shown in fig.[5j 
where the theoretical curves are based on atomic trajectories in the atom-surface 
potential that is assumed to be either of van der Waals-London or of Casimir-Polder 
form. The data are clearly consistent with the CP interaction, hence retardation al- 
ready plays a role for typical distances in the range of 500 nm. 

Subsequent measurements of the Casimir-Polder force have been done, among 
others, using the phenomenon of quantum reflection of ultra-cold atomic beams 
from a solid surface |109[ llO) . In the experiment, a collimated atomic beam of 



metastable Neon atoms impinges on a surface (made of Silicon or some glass) at 
a glancing angle (very small velocity normal to the surface). In this regime, the 
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Fig. 5 (left) Scheme of the experiment of Sukenik & al, taken from Ref. 108 .An atomic beam 
enters a micron- sized gold cavity, and the flux of atoms emerging the cavity is detected and related 
to the atom-surface potential inside the cavity, (right) Measurement of the atom-surface interaction 
in the Casimir-Polder regime, in the experiment of 1 108 1, taken from the same paper. The opacity 
is proportional to the number of atoms which do not exit from the cavity, and is related to the atom- 
surface potential. The solid lines are the theoretical prediction based on: (a) full DLP potential, (b) 
van der Waals-London (short-distance) potential, and (c) no atom-surface potential. 



de Broglie wave of the incident atoms must adapt its wavelength to the distance- 
dependent potential, and fails to do so because the potential changes too rapidly on 
the scale of the atomic wavelength. This failure forces the wave to be reflected, a 
quantum effect that would not occur in an otherwise attractive potential. In the limit 
of zero normal velocity (infinite wavelength), the reflection probability must reach 
100%. The variation with velocity depends on the shape of the atom-surface poten- 
tial and reveals retardation effects | |121[|122| | In fi g. [6|is shown the measurement 
of quantum reflection performed by Shimizu |l09p . In this case, the accuracy was 
not high enough to distinguish reliably between theoretical predictions. More recent 
data are shown in Fig. 5 of the chapter by DeKieviet et al. in this volume. 



Fig. 6 Reflectivity as a function of the nor- 
mal incident velocity of Ne* atoms on a 
Si( 1,0,0) surface, taken from Fig.3 of (109). 
The experimental points (squares) are plot- 
ted together with a theoretical line calculated 
using the approximate expression Vcp = 
-C^l{[d^a)d\ where C4 = 6.8 lO'^^ 
and <3 is a fitting parameter. 
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^ Reprinted figures with permission from C. I. Sukenik, M. G. Boshier, D. Cho, V. Sandoghdar, and 
E. A. Hinds, Phys. Rev. Lett. 70, 560 (1993) Copyright (1993) by the American Physical Society. 
^ Reprinted figure with permission from F. Shimizu, Phys. Rev. Lett. 86, 987 (2001) Copyright 
(2001) by the American Physical Society. 
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The crossover between the van der Waals and Casimir-Polder regimes has been 
recently measured by the group of J. Fortagh |123| , using the reflection of a cloud 
of ultracold atoms at an evanescent wave atomic mirror. This experiment improves 
previous data obtained by the A. Aspect group 1 117] into the crossover region. The 
data are shown in Fig |7] where "vdW" and "ret" label the asymptotes van der Waals 
and Casimir-Polder potentials, respectivel}j^ The full calculation (DLP theory) is 
labelled "trans" and shows some deviation in the crossover region. The data (shown 
with error bars) are clearly favoring the theory including retardation. 




distance (nm) 



Fig. 7 Measured and theoretical prediction 
for the Casimir- Polder interaction, taken 
from fig. 3 of |124|. In the large figure 
theoretical calculation: asymptotic van der 
Waals-London (vdW), asymptotic Casimir- 
Polder (ret), and full theoretical curve 
(trans). In the inset, measured data points 
are included: statistical and systematic errors 
are indicated by the error bars, and the gray 
shaded area, respectively. 



5.3 The Cornell experiments 
5.3.1 Lifshitz regime 

The atom- surface interaction in the Lifshitz regime has been explored in Cornell's 
group 1 5 125 1 . Here a quantum degenerate gas in the Bose-Einstein condensed phase 
has been used as local sensor to measure the atom- surface interaction, similar to the 
work in V. Vuletic's group where smaller distances were involved ||111| . The Cornell 
experiments use a Bose-Einstein condensate (BEC) of a few 10^ ^^Rb atoms that are 
harmonically trapped at a frequency COtrap- The trap is moved towards the surface of 
a sapphire substrate, as illustrated in Fig [8] 

Center-of-mass (dipole) oscillations of the trapped gas are then excited in the 
direction normal to the surface. In absence of atom-surface interaction, the fre- 
quency of the center-of-mass oscillation would correspond to the frequency of 
the trap: (Ocm = cOtmp- Close to the substrate, the atom- surface potential changes 
the effective trap frequency, shifting the center-of-mass frequency by a quantity 



^ Reprinted figure with permission from H. Bender, P. W. Courteille, C. Marzok, C. Zimmermann, 
and S. Slama, Ph ys. Rev. Lett. 104, 083201 (2010) Copyright (2010) by the American Physical 
Society. 
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Fig. 8 (left) Scheme illustrating the experimental configuration in the E. Cornell group, taken 
from Fig.l of 1 125 |. (a) Typical arrangement of the condensate close to the surface. The cloud is 
formed by a few hundred thousand Rubidium 87 atoms, its axial length is ~ 100/im. The surface 
is made of fused silica. The coordinate axes and the direction of gravity are indicated, (b) Typical 
data showing the center-of-mass (dipole) oscillation (x-direction normal to the surface). This is 
obtained after holding the BEC near the surface and then shifting it rapidly away from the surface; 
the "expanded position" is proportional to the velocity component V;^- 

(right) Measured and theoretical frequency of the BEC center-of-mass motion, relative to the nom- 
inal trap frequency ft)trap and normalized as 7 = {(0cm — Ohmp) / ft^trap. Each data point represents a 
single measurement, with both statistical and systematic errors. The mean oscillation amplitude 
is ~ 2.06/im, and the typical size of the BEC (Thomas-Fermi radius) in the oscillation direction 
is ^ 2.40jUm. Theory lines, calculated using theory from 1 126], consider the full atom- surface 
potential: 7 = K (dashed, black line), T = 300 K (solid, blue line) and 7 = 600 K (dotted, red 
line). No adjustable parameters have been used. The result of the van der Waals-London potential 
has been added (dash-dotted, green line). 



7 = {o)cm — 0)trap) / o^rap- The value of J Is related to the atom- surface force |T| and 
for small oscillation amplitudes we have: 

coL = «trap + ^ r dz nliz) dj.^{z) (89) 

III' J — 

where m is the atomic mass, ^{z) is the atom- surface free energy in Eq. ( [2T] ) (z 
is the direction normal to the surface) , and (z) is the xj-integrated unperturbed 
atom density profile Q that takes into account the finite size of the gas cloud. In the 
Thomas-Fermi approximation 

*> = Ti:i'-(i)^ « 

where (typically of few microns) is the cloud radius along z, which depends on 
the chemical potential. In the comparison with the experiment, non-linear effects 
due to large oscillation amplitudes ||4J may become relevant | |125| . The experiment 



Fluctuation-induced forces between atoms and surfaces: the Casimir-Polder interaction 



37 



of Ref . | |125| | was performed at room temperature and succeeded in measuring the 
atom-surface interaction for the first time up to distances L ^ 7/im (see FigjSj). Al- 
though the relative frequency shift in Eq.(89) is only ~ 10~^, the damping of this 
dipole oscillator is so weak that its phase can be measured even after hundreds of 
periods, see FigjSjleft). The same technique has been recently proposed to test the 
interaction between an atom and a non-planar surface 1 127 128| . 



5.3.2 Temperature dependence and non-equilibrium force 



The experiment of Ref. |125| did not reach the accuracy to discriminate between 
the theoretical predictions at T = OK and the T = 300 K, and a clear evidence of 
thermal effects was still missing. In this experiment there was no room to increase 
the temperature of the surface: at high temperatures atoms thermally desorb from 
the walls of the cell, the vacuum in the cell degrades, resulting in the impossibil- 
ity to produce a BEC. To overcome this experimental limitation a new configura- 
tion was studied, where only the surface temperature was increased: the quality of 
the vacuum was not affected because of the relatively small size of the substrate. 
The non-equilibrium theory of atom-surface interactions in this system was devel- 
oped in Refs. | [94l[97l|105||106||129| , as outlined in Sec|4] It predicts new qualita- 



tive and quantitative effects with respect to global equilibrium that are illustrated 
in Fig. fright). The experimental measurement has been achieved in 2007 1 5 1 and 
remains up to now the only one that has detected thermal effects of the electromag- 
netic dispersion interactions in this range of distances. A sketch of the experimental 



apparatus is given in Fig.p] the experimental results in Fig. 10 




Fig. 9 Scheme of the experiment of Ref. j5j 
Pyrex (from which the figure is taken), where atom- 
Holder surface interactions out of thermal equilibrium 

have been measured. 



5.3.3 Outlook 

Precision measurements of the atom-surface interaction may shed light on the on- 
going discussion about the temperature dependence of dispersion interactions with 
media that show absorption, like any conducting medium. It has been pointed out 
by Klimchitskaya and Mostepanenko | |130| that if the small, but nonzero conduc- 
tivity of the glass surface in the Cornell experiment |5] |125| had been taken into ac- 
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Fig. 10 Measured and theoretical frequency 
shift 7 of the center-of- mass frequency cOcm 
for a trapped atomic BEC (^^Rb atoms) close 
to surface (fused silica), in a system in and 
out of thermal equilibrium. The figure is 
taken from fig.4 of |5 1. (a) The figure shows 
three sets of data and accompanying the- 
oretical curves with no adjustable parame- 
ters for various substrate temperatures. Data 
are shown for different substrate tempera- 
tures: Ts = 310 K (blue squares), Ts = 479 
K (green circles), and Ts = 605 K (red trian- 
gles). The environment temperature is main- 
tained at 7e = 310 K. The error bars rep- 
resent the total uncertainty (statistical and 
systematic) of the measurement, (b) Average 
values 7 over the trap center- surface separa- 
tions of 7.0, 7.5, and 8.0 /im, plotted versus 
substrate temperature. It is evident a clear in- 
crease in strength of the atom- surface inter- 
action for elevated temperatures. The solid 
theory curve represents the non-equilibrium 
effect, while the dash-dotted theory curve 
represents the case of equal temperatures. 
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count, the Lifshitz-Keesom tail would involve an infinite static dielectric function, 
and hence deviate from a dielectric medium where 1 < £{0) < oo. This theoretical 
prediction would also be inconsistent with the data. This issue is related to similar 



problems that arise in the macroscopic Casimir interaction, see Refs. |131f[133| for 



reviews. In the atom- surface case, Pitaevskii has pointed out that a smooth crossover 
from a metal to a dielectric is obtained within a non-local description that takes into 
account electric screening in the surface (wave-vector-dependent dielectric function 
e(k, CO)) 1 134] which Ref. |130J did not include. 

New interesting experimental proposals have been presented in order to measure 
the atom-surface force with higher accuracy, essentially based on interferometric 
techniques. All of them deal with atoms trapped in a periodic lattice made by laser 



beams ("optical lattice" 1 135 136|), and placed close to a substrate. Gradients in 



the potential across the lattice can be detected with coherent superposition states of 
atoms over adjacent lattice sites |7 |. These gradients also induce Bloch oscillations 
through the reciprocal space of the lattice: if hq is the width of the Brillouin zone, 
the period Tb of the Bloch oscillations is 

Tb nq 
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where the average (overbar) is over the cloud size in the lattice. The atom-surface 
interaction would, in fact, only shift the Bloch period by a relative amount of 
10~^ . . . 10~^ if the main force is gravity and the atoms are at a distance L ^ 5 /im 
| |1Q6| . Distance-dependent shifts in atomic clock frequencies have also been pro- 
posed |137| . They arise from the differential energy shift of the two atomic states 
which are related to the difference in polarizabilities. Finally, a corrugated surface 
produces a periodic Casimir-Polder potential that manifests itself by a band gap 
in the dispersion relation of the elementary excitations of the BEC (Bogoliubov 
modes). The spectrum of these modes is characterized by a dynamic structure factor 
that can be detected by the Bragg scattering of a pair of laser beams |138| . 
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